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We propose to use a time dependent imaginary potential to describe quantum mechanical tunnelling 
through time- varying potential barriers. We use Gamow solutions for stationary tunneling problems 
to justify our choice of potential, and we apply our method to describe tunneling of a mesoscopic 
quantum variable: the phase change across a Josephson Junction. The Josephson Junction phase 
variable behaves as the position coordinate of a particle moving in a tilted washboard potential, and 
our general solution to the motion in such a potential with a time dependent tilt reproduces a number 
of features associated with voltage switching in Josephson Junctions. Apart from applications as 
artificial atoms in quantum information studies, the Josephson Junction may serve as an electric field 
sensitive detector, and our studies provide a detailed understanding of how the voltage switching 
dynamics couples to the electric field amplitude. 
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I. INTRODUCTION 



Tunneling is a widely observed phenomenon in quantum 
mechanics and for tunneling through stationary barriers, 
scattering theory and good approximations based on the 
JWKB method are available pQ. For the general dynam- 
ical case such simple approximations do not apply and 
only in certain cases is it possible evaluate the tunnel- 
ing rate [2]. The part of the wave function which is still 
trapped behind the barrier is evolving within the trapping 
region, and it may repeatedly encounter the tunnel barrier 
in a highly non-trivial manner. While the Schrodinger 
equation is uniquely defined in a time dependent potential, 
the handling of the scattering continuum components of 
the wave function is impeded by the demand for precise 
calculations on a very large interval of the coordinate 
variable. In this paper we propose a new, effective ansatz 
to handle tunneling in a time dependent potential. The 
method addresses tunneling in a general manner, but as 
a key example we will focus our attention on tunneling 
of a mesoscopic variable, the phase of a current-biased 
Josephson junction(CBJJ). 

Circuit quantum electrodynamics components have re- 
cently reached an experimental state of perfection that 
makes quantum interaction and control properties of these 
man-made devices compete with, or supersede the ones of 
natural, microscopic quantum systems, while at the same 
time retaining their fundamental interest as macroscopic 
quantum systems [3l [4] . The CB J J is a device of par- 
ticular interest since the phase variable associated with 
the superconductor is trapped in a washboard potential, 
and its quantum level structure permits identification and 
control of pairs of states, which may be used as qubits for 
the realization of quantum information processing [5] . As 
shown in Fig.l, its motional states are subject to tunnel- 
ing which corresponds to a voltage drop over the device 
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FIG. 1. (Color online) The superconducting phase variable 
in a current-biased Josephson junction behaves as a quantum 
particle trapped in a tilted washboard potential. The dashed 
(green) line is a sketch of the ground state wavefunction in 
one of the trapping wells. The probability that the phase 
particle tunnels out of the potential is increased if the barrier is 
modulated, e.g., by variation of the tilt or by a time dependent 
microwave field. 

and which may be used to distinguish the qubit states. 
Also, by modulation of the potential tilt or by a radiation 
field, a given initial state may be driven into states that 
undergo tunneling to running states. The switching can 
occur due to thermal activation [6] and due to macroscopic 
quantum tunneling. In this treatment the contribution 
from thermal activation will be neglected. This implies 
that the system can be described as an ordinary quantum 
particle and general methods from quantum mechanics 
can be applied to qualitatively and quantitatively study 
several features associated with its dynamics. While our 
approach to tunneling in a time dependent potential is 
general, we will here discuss its application to the CBJJ 
and present results analogous to observations made in 
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CBJJ experiments. 

The CBJJ has also been suggested as a sensitive mi- 
crowave photodetector [8] , and as a first step towards a 
low frequency quantum detection theory for such a device, 
we will study its response to driving by a weak microwave 
field. 

A weak field-induced modulation, resonantly coupling 
the initial state to eigenstates subject to stronger tunnel- 
ing, can be treated as a perturbation and for the CBJJ 
theory yields good agreement with experimental results 
[9HT2]. In this paper we will describe the wave- function 
subject to a strongly modified potential. This case is not 
amenable to the perturbative description, and instead, 
we shall formally describe the tunneling process by wave 
packet propagation in a time dependent imaginary poten- 
tial (TDIP), such that tunneling, observed as switching of 
the CBJJ into the voltage state, is treated as a dynamical 
loss process. In this way we both determine the time 
dependent probability for the tunneling to occur and we 
appropriately describe the evolution of the wave function 
conditioned on no tunneling event being detected. 

The paper is organized as follows. In Sec. [IT] we define 
the tunneling problem, we motivate our TDIP description, 
and we provide an explicit time dependent expression for 
the imaginary potential applied in our numerical studies. 
In Sec. |III| we use our method to study different regimes of 
tunneling dynamics with increasing bias currents applied 
to the CBJJ. In Sec. [lV]we fix the bias current and we 
analyze the performance of the Josephson junction as a 
field detector by examination of the switching properties 
under application of a driving field. Sec. [V] concludes the 
paper. 



II. 



TUNNELING OF THE JOSEPHSON 
JUNCTION PHASE 

Tunneling loss and complex absorbing 
potentials 



The Josephson junction phase (j) (see Sec II B) behaves 
as the position of a particle described by the time depen- 
dent Schrodinger equation, 



d 



(1) 



where the Hamiltonian operator H(t) contains a poten- 
tial and an effective kinetic energy term. The tunneling 
dynamics of the phase variable, which physically corre- 
sponds to the switching between different voltages across 
the device, is hence mathematically described as a con- 
ventional quantum tunneling problem. The purpose of 
this manuscript is to establish an effective, approximate 
theory for the tunneling dynamics in a time dependent 
potential, with special attention to parameter regimes 
relevant for the CBJJ. Following the early description by 
Gamow [1], one can describe tunneling in static poten- 
tials by eigenfunctions of the time-independent Hamilton 



operator. If these functions are chosen with outgoing, ra- 
diating boundary conditions and used as a wave function 
basis, we effectively redefine the scalar product on the 
system Hilbert space, and the original Hamiltonian is not 
Hermitian. This in turn leads to the emergence of com- 
plex eigenvalues, and the loss of norm associated with the 
imaginary part of the energy eigenvalues represents the 
probability of tunnelling [13 -15 . Effective non-Hermitian 
Hamiltonians also emerge in quantum optics and quantum 
measurement theory, where they govern the evolution of a 
quantum system conditioned on the absence of absorption 
or loss events. In these theories, the decreasing norm of 
the state vector provides the probability for the evolution 
to occur without these events happening, while one may 
simulate the complete dynamics including random detec- 
tion events by suitable application of "quantum jump" 
operators [T6HT8] . 

At a given time the tunneling probability is governed by 
the wave function obtained by propagation with our non- 
Hermitian Hamiltonian, i.e., the state conditioned on the 
detection of no previous switching events. The amount of 
trapped population and its actual wave function therefore 
yields the local loss rate due to tunneling and the whole 
time dependent dynamics should be well described by the 
Gamow vectors and their complex eigenvalues. Another 
approach to yield complex eigenvalues is by directly imple- 
menting a non-Hermitian Hamiltonian with respect to the 
original Hilbert space of square integrable wave functions, 
for example by the introduction of a complex absorbing 
potential (CAP). A properly designed CAP may thus 
describe the same essential physics as the Gamov vectors, 
i.e., the evolution of the unnormalized states may yield 
identical or very similar wave function behaviour in the 
spatial and temporal range of interest [19] . 

There is a rich literature [20] on the identification of 
suitable CAPs, but since we shall be dealing with the 
further complexity of tunneling through a time dependent 
potential, we shall merely propose a simple, physically 
motivated ansatz for our time dependent imaginary po- 
tential (TDIP), Vi m ((j),t), and solve the time-dependent 
Schrodinger equation 

ih^(4>, t) = H (i)V(^, t) - iv im (<p, t)ip(<p, t). (2) 



B. The dynamics of a CBJJ 

In this section we outline the particular properties and 
the effective quantum description of a CBJJ by a one- 
dimensional Schrodinger equation for the phase variable 
(j). The Josephson junction is composed of two supercon- 
ductors, R(ight) and L(eft), separated by a thin isolating 
layer. The macroscopic wave functions of R and L, i/jr^l 
differ by a phase factor ipL = e^ipR, which constitutes 
the macroscopic quantum degree of freedom of the system. 
This phase can be described [21] as a quantum particle 
with mass M = C(^>o/27r) 2 , where C is the capacitance 
of the junction and &o = 7^ is the flux quantum. The 
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particle experiences a potential 



U (<j>) = -Ej(I<j> + cos 4>) 



(3) 



where Ej = is the Josephson energy and / = ^, 

with ^ being the bias current applied to the junction. 
When It exceeds the critical current 7 C , the potential tilt 
dominates the harmonic variation with 0, and the phase 
becomes classically unbound. 

We further include the interaction of the CBJJ with a 
time dependent microwave field via the potential term, 



Urr 



-E 



j770sin (u mw t) 



(4) 



where 77 is proportional to the field amplitude. 

The CBJJ is thus described by the time-dependent 
Schrodinger equation 



C/(0,t)W(0,t) (5) 
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with U((j),t) = Uo + Umw- This time dependent potential 
is sketched in Fig. 1, indicating also the tunneling process, 
responsible for the switching of the Josephson junction. 



C. Outgoing states and absorbing potential 

An imaginary absorbing potential, extending beyond 
the outer turning point of a potential barrier would seem a 
natural candidate to remove the tail of the wave function 
as tunneling develops. Since we want the potential to 
remove the projection of our wave function on the running 
states in that region, it is useful to apply approximate 
solutions for their position (phase) and time dependence 
in our ansatz for the TDIP. From [22] we have an approx- 
imate expression for the running state of the Josephson 
junction phase variable > (j)turm initially at time r in 
the ground state, in the linear potential region beyond 
the outer classical turning point of the potential barrier 



Our aim is to remove the projection of our time de- 
pendent wave function \ij){t)) on a suitable set of running 
states of the form P ou M) = Z)^ JVw)(Vw|^) • 
We model this operation by a time-dependent complex 
potential, obtained as an integral over different running 
states, emitted within the past interval of time At, 



\^out(M)\ 2 dt' . 



(7) 



t-At 



The parameter (3 serves both as an adjustable strength 
parameter and as a convenient normalization for the tem- 
poral integration. We assume a constant value of T and 
At, with At -1 <C r <C cop,, and we evaluate the integral 
([7]) only for values larger than the outer turning point, 
<t>tum, of the potential. The dependence on T drops out, 
and we obtain 



for 



> 
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for < <j>turn(t)- 



(8) 



Here, we recall the time-dependence of the plasma fre- 
quency as the effective bias-current is changing with time. 

Unlike normal CAPs used in time-independent prob- 
lems, this potential is modelled to absorb the running 
state components |6| pertaining to the time dependent 
Hamiltonian, and our ansatz imaginary potential indeed 
attains finite values only beyond the time dependent outer 
classical turning point of the real potential. We emphasize 
that ([8| is only an ansatz, but it ensures that the deple- 
tion of the wave function norm is properly associated with 
the probability that the particle has tunnelled through 
the barrier. 

While it is generally a challenge to design CAPs that 
do not reflect part of the wave packets impinging on them, 
our application only concerns the small wave function 
amplitude that has tunneled through the real potential 
barrier, and we have verified that reflection is insignificant 
in our calculations. To account for the running solution, 
the TDIP must furthermore suppress the wave function 
before it hits and gets reflected by the border of the grid 
used for the calculation, and a sufficiently large grid is 
readily identified in the numerical calculations. 
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with 0' = — 0o, where 0o is the initial equilibrium posi- 
tion (bottom of the well) and hu is the energy difference 
between the bottom of the well and the energy of the 
lowest quasi-bound state. We have further introduced 

the frequency parameter ujq = yj 2 j^ and the plasma 

frequency uo p = ujq{1 — I 2 ) 1 / 4 . The wavefunction ^ is 
defined for tuo p > \/20; and we assume / < 1. The rate 
parameter T in (|6| which attains a definite value in a 
static potential, will be briefly discussed below. 



D. Friction and junction resistance 

We have now presented our candidate TDIP, but be- 
fore we proceed to numerical examples let us include the 
effects of friction (junction resistance). Methods have 
been developed to incorporate friction and diffusion in 
the quantum mechanical description of Brownian mo- 
tion of a particle coupled to its environment, [23l [24] . 
Assuming zero temperature and Markovian noise corre- 
lation in the environment, we can describe the effects 
of the junction resistance, R, by a non-linear imaginary 

potential term, —ic(4> ~ (fyt) -> where ( ex and 
(0)t = |'0(0, t)\ 2 dej) is the mean value of at time 
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FIG. 2. (Color online) Blue (solid) bars show the switching 
current distribution simulated with an absorbing potential. 
Red (hollow) bars are calculated by Eq. (10). We have chosen 
C and R such that uj q = 0.0183 Ej/h and C = 8.4 x 10~ 4 £j. 
The bias current is increased linearly from / = 0.2 at a slow 
rate, dl/dt = Results are shown for / > 0.95. 



t. When added to the Hamiltonian this term penalizes 
large variations of <p around its mean and decoheres the 
spatial wave function, which in turn leads to friction for 
the phase variable <f>. It also leads to a loss of norm, albeit 
on a typically slower scale than the tunneling dynamics. 
We shall include this term in our calculations on the CBJJ, 
but we shall renormalize the wave function with respect to 
the loss it incurs, so that we unambiguously associate the 
wave function loss of norm with the tunneling dynamics. 

To summarize, we treat the entire problem by solving 
the non-linear Schrodinger equation 



d 
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iV im ((j),t)ilj((j),t) 



<(0-(0>t)%(0,t) (9) 



and renormalize with respect to the loss of norm caused 
by the last term. 



E. Tunneling rates and and switching current 
distributions 



In a time-independent potential, the tunneling rate of 
the Josephson junction phase variable has been deter- 
mined by Caldeira and Leggett [25) 126], 
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where AU = 2Ej(y/l — I 2 — 7acos(/)) is the barrier- 
height, and where effects due to the friction are also taken 



explicitly into account. This analytic expression available 
for CBJJs qualify them as ideal candidates to test our 
general tunneling description. 

The gradually decreasing norm, HV'WII 2 of our nu- 
merically determined wave packet is interpreted as the 
probability that the phase variable has not tunneled until 
time t. The probability for a current switching event in 
the next infinitesimal time interval dt, is thus simply given 
by the loss of norm in that interval, and conditioned on 
no previous event, the switching rate of the CBJJ reads 



It 



d\\m\? 

dt 



iwmw 2 . 



(ii) 



It is convenient in experiments to determine the switch- 
ing current distribution, i.e., the probability distribution 
P(I) for switching events to occur at different values of 
the bias current /, while this is being ramped up slowly, 

i(t) = i + f t .t. 

Our wave packet propagation yields the surviving (non- 
switching) population HVK^H 2 as function of time, and 
we directly obtain the corresponding switching current 
distribution, obeying 



P(I)AI 



dt 

d\mw 

dt 



At 



■dl\- 
.dt) 
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evaluated at the time t such that I = I(t) and At being 
the infinitesimal time interval in which / increases by AI. 

Under the assumption of a slowly ramped bias current, 
the rate jcl found by Caldeira and Leggett leads to a 
switching current density, expressed as a product of the 
current tunneling rate with the survival probability until 
the value I is reached during the ramp, 



CL 



dJy 
~dt) 



lCL(I)e~ 



(13) 



In Fig. [2] we see, that our calculation, matches the 
result of the quasi-static switching current distribution, 
(13), reasonably well as a function of the bias current. 



III. DRIVING WITH INCREASING BIAS 
CURRENT 

To study a system with a more complicated time- 
dependent tunneling dynamics we now include a mi- 
crowave field with a constant power (constant 77) and 
a frequency uu mw , while we increase / with a low constant 
rate. Many experimental results are available for the 
CBJJ under such conditions, and this special example of 
tunneling dynamics serves as the primary application of 
our theoretical description. 

The evolution of the system is calculated numerically 
using equation (|9| with the imaginary potential given in 
(pi). The parameters are chosen to represent a Josephson 
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FIG. 3. (Color online) Switching current distributions for a CBJJ driven by a field. In all figures we have chosen C and R such 
that cjo = 0.0183 Ej/h and £ = 8.4 x 10~ 4 Ej. The frequencies and hence the resonance currents are chosen as u) m w = 0.64cjo, 
I — 0.90 in (al) and (a2), uj mw — 0.51cjo, I — 0.95 in (bl) and (b2), uJmw — 0.36cjo, I — 0.98 in (cl) and (c2), and uj mw — 0.255o;o- 
/ = 0.987 in (dl) and (d2). In panels (cl-2) and (dl-2) the resonance current is above the current of the primary zero-field peak. 
In (al-dl) the applied microwave amplitude attains the values, r\ = 0, 0.002, 0.004, 0.006, 0.008, 0.010, 0.014, 0.018, 0.022, 0.030, 
0.040, 0.055, 0.070, 0.090, 0.115, 0.140, which yield the curves shown in the front (dark-blue) towards the back (dark-red) in 
each panel. In (a2-d2) we see the color plots showing the same result as the three dimensional plots in (al-dl). The calculations 
are carried out with a bias current that increases linearly with time from / = 0.2 at a rate equal to dl/dt = (In (al-dl) 

the results are only shown for / > 0.6, since P(I) <C 1 for the smaller bias currents. In (a2) and (c2) we show the results for 
/ > 0.7, while in (b2) and (d2) the results are shown for / > 0.8.) 



junction with a critical current of ~2 A. The results of 
the simulations are shown in Fig. [3^al-dl) where the 
probability distribution is plotted as three dimensional 
plot, while we in [3^a2-d2) plot the switching current 
probability distributions as two dimensional color plots to 
give both a qualitative and a quantitative overview of the 
results. In the following subsections we will separately 
discuss many of the features observed in Fig. 3. 



A. Resonance Peaks 

We notice in Fig. [3jal-2) and[3jbl-2) that when the 
microwave field strength is increased, in addition to the 
conventional switching current peak a second peak ap- 
pears. In the two dimensional color plots this is seen 
as the opening of a gap in the switching current curve. 
This second peak is associated with the tunnelling via 
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C. Dynamical Bifurcation 
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FIG. 4. (Color online) Switching current distributions for a 
CBJJ driven by a field. The parameters are chosen as in Fig. 
[3] (a) and with 77 = 0.008. These parameters give rise to a 
multi-peak structure of the switching current distribution. 



the excited phase eigenstate [27j [28] , and hence the peak 
starts appearing where the frequency uo mw = ojqi(I), with 
the resonance frequency 



(14) 



In Fig. [3^dl-2) at a much higher field strength we 
observe a second peak at the same bias current as in 
Fig. [3jbl-2). But, since, in this panel, the microwave 
frequency is half the resonance frequency, the excitation 
of the phase variable is a second order process [29] . 

When the power is further increased the peaks move lin- 
early to lower bias current, which is also well-understood 
[30] . This behaviour is clearly observed in Fig. [3|a2-d2). 
The strong microwave field then effectively suppresses the 
potential barrier as the power is increased. 



B. Multi-peak Structure 

In Fig. |3ja), after the revival of the second peak, a 
multi-peak structure appears. This is shown in Fig. [4] 
Since we only drive the junction with a single frequency 
field, these extra peaks cannot be understood as con- 
ventional resonances. They are, however, a well known 
effect associated with the multi-peaked Fourier transform 
of a frequency chirped electric field amplitude [31], see 
also [32]. The analysis in [31 leads to peaks with higher 
density for the largest detunings and, although the cor- 
respondence between the two physical problems is only 
approximate, it is consistent with our observation that 
the multi-peak structure is not discernible in Fig. |3|b). 
We should emphasize that in typical experiments, ^ is 
smaller than in our calculations leading to more narrow 
peaks than seen here. 



In Fig. [3^cl-2) we see yet another interesting feature 
of a driven CBJJ. The resonance condition lies here at 
a bias current above the zero-field peak current, which 
implies that we see no resonance peak. However for very 
strong fields we observe a very broad splitting into a multi 
peak structure. Such results are usually explained by a 
dynamical bifurcation [27] [33] . Further analysis of this 
phenomenon is beyond the scope of the present work, but 
it is reassuring that our simple model also gives rise to 
this complex behaviour. 



IV. DRIVING WITH CONSTANT BIAS 
CURRENT 

We now maintain the bias current at a constant, 
/ = 0.96, close to the switching value, and we apply 
a microwave field to the junction. As we have already 
seen, the tunnelling rate is increased if the resonance 
condition is met. In Fig. [5] we observe what happens 
at different frequencies. First, the switching probability 
follows the field strength, but as we approach resonance, 
a peak emerges, and when the frequency is increased even 
further the peak is reduced. In the color plot we see a 
hyperbolic behaviour of the signal following each oscilla- 
tion cycle, and for increasing frequency they, naturally, 
move closer together and increase in strength and as a 
consequence the later ones disappear. After the resonance 
we see that the signal slowly smears out. 

A switching rate that follows the field strength in real 
time allows implementation of a field amplitude detector 
capable of resolving single oscillations of the field. In 
Fig. [6] we increase the field strength while maintaining a 
very low frequency and see that when the field is strong 
enough single positive oscillations of the field amplitude 
are clearly detected. In the color plot we also clearly 
see the growth and reduction of equally spaced areas of 
probability. 

We intend to study the performance and sensitivity of 
this detection mode of the CBJJ further with particu- 
lar emphasis on the prospects for quantum information 
processing in circuit quantum electrodynamics. 



V. CONCLUSION AND OUTLOOK 

We have introduced a novel, effective method using time 
dependent imaginary potentials to describe the switching 
behaviour of a current biased Josephson junction. The 
method is readily implemented with standard wave packet 
solvers and with realistic parameters, it reproduces a wide 
range of results that have previously been investigated 
experimentally and theoretically by other techniques [27]- 
[33] . We emphasize that our treatment builds on an ansatz 
for a time dependent imaginary potential (TDIP), and a 
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FIG. 5. (Color online) Switching distributions for different frequencies with a fixed bias current at / = 0.96. We have chosen 
C and R such that u = 0.0183 Ej/h and C = 8.4 x 10~ 4 £j. This gives a zero-field energy splitting between the ground and 
first exited states at cjoi = 0.48cjo- To the left, we have increasing frequency from the front (dark-blue) to the back (dark-red) 
at ujmw = (0.07, 0.23, 0.38, 0.54, 0.69, 0.84, 1.00, 1.15, 1.31)cj i. To the right, we show a color plot of the same switching 
distributions where the signal from later cycles decrease and eventually disapper. The power of the field is chosen with r\ = 0.020. 
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FIG. 6. (Color online) Switching distributions for different amplitudes with a fixed bias current at I — 0.96. We have chosen 
C and R such that cjo = 0.0183 Ej/h, cjoi = 0.48cjo and f = 8.4 x 10 -4 i^j. We drive the junction at a very low frequency 
oJmw = 0.036o;o. To the left the power of the field is increasing from the front (dark-blue) to the back (dark-red) at 77 = 0.012, 
0.016, 0.020, 0.024, 0.028 and 0.032. The arrows (magenta) indicate when the accumulated probability for a switch equals 
unity. In the strong field case this happens after only one half cycle. To the right we show a color plot of the same switching 
distributions as a function of the driving field the amplitude. 



number of possibilities may be explored for quantitative 
improvement of the potential chosen. 

An advantage of the TDIP method is its ability to 
deal with explicitly time dependent driving fields. With 
the explicit treatment of the wave function of the phase 
variable conditioned upon the tunneling dynamics, we 
may readily extend the studies towards other systems 
where non-trivial tunneling dynamics in time-dependent 
potentials is present. We may also expect that the meth- 
ods developed here can be used to model the dynamics 
of coupled CBJJs and of correlated and entangled CBJJ 



dynamics with the quantized radiation fields. Potentially 
our approach may thus form the basis for a novel quantum 
theory of measurement, where the noise and back action 
is adapted to the special field amplitude sensitivity of the 
device. 
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